Kondo effect in the helical edge liquid of the quantum spin Hall state 
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Following the recent observation of the quantum spin Hall (QSH) effect in HgTe quantum wells, 
an important issue is to understand the effect of impurities on transport in the QSH regime. Using 
linear response and renormalization group methods, we calculate the edge conductance of a QSH 
insulator as a function of temperature in the presence of a magnetic impurity. At high temperatures, 
Kondo and/or two-particle scattering give rise to a logarithmic temperature dependence. At low 
temperatures, for weak Coulomb interactions in the edge liquid the conductance is restored to uni- 
tarity with unusual power-laws characteristic of a 'local helical liquid', while for strong interactions 
transport proceeds by weak tunneling through the impurity where only half an electron charge is 
transferred in each tunneling event. 

PACS numbers: 71.55.-i, 72.25.Hg, 73.43.-f, 75.30.Hx 
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The quantum spin Hall (QSH) insulator is a topologi- 
cally non-trivial state of matter |l| that has recently been 
observed in transport experiments carried out in HgTe 
quantum wells (QW) following its theoretical predic- 
tion Q . The two-dimensional QSH insulator has a charge 
excitation gap in the bulk but supports one-dimensional 
gapless edge states forming a so-called helical liquid (HL) : 
on each edge there exists a Kramers' pair of counter- 
propagating states with opposite spin polarization. The 
QSH insulator is robust against weak single-particle per- 
turbations which preserve time-reversal symmetry such 
as weak potential scattering 0, 0, Q . 

This theoretical picture is consistent with experimen- 
tal observations: the longitudinal conductance G in a 
Hall bar measurement is approximately quantized to 
Go = 2e 2 /h, independent of temperature, for samples 
of about a micron length @, 0|- However, larger samples 
exhibit G < Gq and G decreases with decreasing temper- 
ature Q ■ Deviations from the expected quantized value 
have been attributed to the presence of local doped re- 
gions due to potential inhomogeneities within the sample 
arising from impurities or roughness of the well/barrier 
interface Q. Although pure potential scattering cannot 
backscatter the edge states, the role of these potential in- 
homogeneities is to trap bulk electrons which may then 
interact with the edge electrons. These localized regions 
act as dephasing centers for the edge channels due to 
interaction effects and may cause backscattering. 

In this work, we study theoretically the temperature 
dependence of the edge conductance of a QSH insulator. 
We consider the case where a local doped region in the 
vicinity of the edge contains an odd number of electrons 
and acts as a magnetic impurity coupled to the helical 
edge liquid. Our main results are as follows (Fig. [T|): 

1. At high temperatures, the conductance is logarith- 
mic, -AG ee —{G -Go) = V + lH D / T ) wh ere 77,7 
are interaction-dependent parameters and D is an energy 
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FIG. 1: Temperature dependence of the conductance: the 
behavior is logarithmic at high temperature T ^> T* and 
power-law at low temperature T<T'. At T — 0, a metal- 
insulator quantum phase transition is driven by Coulomb in- 
teractions in the helical liquid: the system is a 'Kondo metal' 
for K > 1/4 and a 'Luttinger liquid insulator' for K < 1/4. 
The Fano factor e* is defined as the ratio between shot noise 
and current, and reflects the charge of the current-carrying 
excitations. 



scale of order the bulk gap. 

2. For weak Coulomb interactions K > 1/4 where K 
is the Luttinger parameter of the edge liquid, the con- 
ductance is restored to unitarity at T = due to the 
formation of a Kondo singlet. This is in stark contrast 
with the Kondo problem in a usual spinful ID liquid 
where the conductance vanishes at T = for all K p < 1 
where K p is the Luttinger parameter in the charge sector 

At low but finite T the conductance decreases as an 
unusual power-law AG cx — T 2 ( 4K_1 ) due to correlated 
two-particle (2P) backscattering. The edge liquid being 
helical, the decrease in conductance is a direct measure 
of the spin- flip rate Q . 

3. For strong Coulomb interactions K < 1/4, 2P 
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backscattering processes are relevant and the system be- 
comes insulating at T = 0. At low but finite T, the 
conductance is restored by tunneling of excitations with 
fractional charge e/2 and we obtain G(T) cx T 2 ( 1 /4^-i)_ 
Model. — We model the impurity by a S — \ local 
spin coupled by exchange interaction to the ID HL with 
Coulomb interactions. The HL having the same num- 
ber of degrees of freedom as a spinless fermion, a sin- 
gle nonchiral boson (j> is sufficient for its description in 
the bosonized language The system is described 

by the Hamiltonian H = Hq + Hk + i?2 where Ho 
is the usual Tomonaga-Luttinger Hamiltonian Ho = 
| / dx [KTl 2 + -^(d x (/)) 2 ] , with K the Luttinger parame- 
ter and v the edge state velocity. The Kondo Hamiltonian 
Hk has the form 



J\\a , 



i2Vw<A(0) . 



h.c.) - ^SJI(O), (1) 



where S± = S x ± iS y and S z are the spin operators for 
the impurity localized at x — 0. a is the lattice constant 
of the underlying 2D lattice and corresponds to the size 
of the impurity (we assume that the impurity occupies a 
single lattice site). £, the penetration length of the he- 
lical edge states into the bulk, acts as a short-distance 
cutoff for the ID continuum theory in the same way that 
the magnetic length, the penetration length of the chiral 
edge states, acts as a short-distance cutoff for the chiral 
Luttinger liquid theory of the quantum Hall edge [Io| . 
In addition to Kondo scattering, 2P backscattering is al- 
lowed by time-reversal symmetry @, @] . In HgTe QW the 
wavevector at which the edge dispersion enters the bulk 
is usually much smaller than 7r/2a such that the uniform 
2P backscattering (umklapp) term requiring Akp = 2ir/a 
can be ignored. The impurity potential can however pro- 
vide a Akp momentum transfer and we must generally 
also consider a local impurity-induced 2P backscattering 
term [IJ H 2 = : cos 40W>(O) : where A 2 is the 2P 
backscattering amplitude. 

Weak coupling regime. — We first consider the weak 
coupling regime where Ju , J z and A2 are small parame- 
ters. The calculation of the conductance proceeds in two 
steps. We first perform an explicit perturbative calcula- 
tion of the conductance to quadratic order in the bare 
couplings Jy and A2 using the Kubo formula. This result 
is then extended to include all leading logarithmic terms 
in the perturbation expansion by means of a weak cou- 
pling renormalization group (RG) analysis of the scale- 
dependent couplings J\\(T) and A2(T) where the scale is 
set by the temperature T. 

The forward scattering term J z can be removed from 
the Hamiltonian by a unitary transformation [l2| of the 
form U = e iAS ^°) with A = -J z a/vK^, at the ex- 
pense of modifying the scaling dimension of the ver- 
tex operator : e l2 V^<t> ■ , The transformed Hamiltonian 
reads UHW = H +H 2 +^ (s_ : e -* 2 v^x0(o) . + hx /j ( 



where x = 1 — vJ z /ZK with v = a/nv the density of 
states of the HL. The scaling dimension of : e a ^ x ^ : is 

k = k x 2 . 

Using the transformed Hamiltonian we obtain the cor- 
rection to the conductance AG(T) = AG K (T) + AG 2 (T) 
to quadratic order in the couplings Jy and A2 using the 
Kubo formula, where AGk is the correction due to spin- 
flip Kondo scattering [l3j ]. 



AGk 

e 2 /h 
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1 S(5+l)(^J,|(T)) 2 , (2) 



where vJ\\(T) = uJ\\(T '/ ' D)^ 1 to 0{J\\) and D = hv/n^ 
is a high-energy cutoff of order the bulk gap. AG2 is the 
correction due to 2P backscattering, 



ag 2 = r(i)r(4g) a 4 / A 2 (T) 

e 2 /h ra + AK) 27r 4 e 4 V D 



(3) 



where A2(T) = A 2 (T/_D) 4 - ftr_1 . There are no crossed 
terms of the form 0(J\\ A2) or 0(J Z \2) for a S = \ impu- 
rity since the 2P backscattering operator flips two spins 
but a S = \ spin can be flipped only once. We however 
expect that such terms would be generated for impurities 
with higher spin. 

These results can be complemented by a RG analysis. 
The RG equation for A2 follows by dimensional analysis 
^ = (l-4fT)A 2 with £ = \n{D/T), so that A 2 is relevant 
for K < 1/4 and irrelevant for K > 1/4. The renormal- 
ized coupling is ^{T) — \ 2 {T / D) 4K ^ 1 , and second order 
renormalized perturbation theory AG2(T) cx — A 2 (T) 2 
simply reproduces the Kubo formula result ([3]). Pertur- 
bation theory fails for T < T* where T* cx (A!]) 1 /^ 4 ^ is 
a scale for the crossover from weak to strong 2P backscat- 
tering with A° the bare 2P backscattering amplitude. 

The one-loop RG equations 0, [3] for the Kondo cou- 
plings J\\,J Z read 



dJ\\ 

~dJ 



{1-K)J» +vJ\\J z 



dJ z 

~dJ 



(4) 



The family of RG trajectories is indexed by a single scal- 
ing invariant c = (i/J||) 2 — {vJ z ) 2 where vJ z = i/J z +l—K, 
which is fixed by the couplings at energy scale D. In 
contrast to the spinful case [8(, the absence of spin-flip 
forward scattering in the HL preserves the stability of the 
ferromagnetic fixed line, as is the case in the usual Kondo 
problem. The renormalized spin-flip amplitude Jn (T) is 
given in terms of the bare parameters J9, J z by 



vJ\\{T) = 



sinh(az/J|?ln(T/T*) 



(5) 



such that AGk is obtained to all orders in perturbation 
theory in the leading-log approximation by substituting 

Eq. © in Eq. ©. In Eq. ©, a = [(J z / 'J,?) 2 - l] 1/2 
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FIG. 2: Strong coupling regime: the Kondo singlet effectively 
removes one site from the system, a) Luttinger liquid with 
K p < 1: a punctured ID lattice is disconnected, b) QSH edge 
liquid with K > 1/4: the edge liquid follows the boundary of 
the deformed 2D lattice, c) Half-charge tunneling for K < 1/4 
by flips of the Ising order parameter m. 



is an anisotropy parameter [l5j | and T£ is the Kondo 
temperature, 



T£ = I>exp - 



1 arcsinh a 



a 



In the isotropic case a = 0, one recovers the usual form 



De 



-l/uJ° 



In the limit l — K^> z/J,?, vJ z , Coulomb 



interactions dominate over Kondo physics and we obtain 

1/(1-K) 



D [ ) t a power-law dependence similar 



to results previously obtained [161 ] for Kondo impurities 
in spinful Luttinger liquids. From the scaling exponent 
we see that T£ corresponds to the scale of the mass gap 
opened in a spinless Luttinger liquid by a point potential 
scatterer of strength v Jy and the corresponding crossover 
is that of weak to strong single-particle backscattering. 

In the high temperature limit max{T 2 *, T£} < T < D, 
both the Kondo and 2P scattering processes contribute 
logarithmically to the suppression of the conductance, 
— AG(T) ~ i] + j\n(D/T) where 77,7 are functions of the 
bare couplings K, J,?, J° and A°. 

Strong coupling regime. — We now investigate the low 
temperature regime below the crossover temperatures 
T < min{T 2 *,T£}. The topological nature of the QSH 
edge state as a 'holographic liquid' living on the bound- 
ary of a 2D system [5| results in a drastic change of the 
low-energy effective theory in the vicinity of the strong 
coupling fixed point as compared to that of a usual ID 
quantum wire. As suggested by the perturbative RG 
analysis, the nature of the T = fixed point depends on 
whether K is greater or lesser than 1/4. 

We first consider the case where K > 1/4. In this case, 



2P backscattering is irrelevant and AG2 flows to zero. 
On the other hand, for antiferromagnetic J z the Kondo 
strong coupling fixed point J11 , J z — > +00 is reached at 
T = 0, with formation of a local Kramers singlet and 
complete screening of the impurity spin by the edge elec- 
trons. As a result, the formation of the Kondo singlet 
effectively removes the impurity site from the underly- 
ing 2D lattice (Fig. [2Jd). In a strictly ID spinful liquid, 
this has the effect of cutting the system into two discon- 
nected semi-infinite ID liquids (Fig. [2^) and transport is 
blocked at T = for all K p < 1 [8]. In contrast, due to 
its topological nature, the QSH edge state simply follows 
the new shape of the edge and we expect the unitarity 
limit G — 2e 2 /h to be restored at T = 0. For finite 
T <C ?kj the effective low-energy Hamiltonian contains 
the leading irrelevant operators in the vicinity of the fixed 
point. In the case of spinful conduction electrons, the 
lowest-dimensional operator causing a reduction of the 
conductance is single-particle backscattering. However, 
the helical property of the QSH edge states forbids such a 
term and it is natural to conjecture that the leading irrel- 
evant operator must be the 2P backscattering operator 
with scaling dimension AK. We thus expect a correc- 
tion to the conductance at low temperatures T <^T^ for 
K > 1/4 of the form AG cx - {T/T£) 2(iK ~ 1} . 

In particular, in the noninteracting case K = 1 we 
predict a T 6 dependence in marked contrast to both the 
usual Fermi liquid [17( and spinful ID liquid [SJ behaviors. 
This dependence characteristic of a 'local helical Fermi 
liquid' can be understood from a simple phase space argu- 
ment. The Pauli principle requires the 2P backscattering 
operator to be defined through a point-splitting proce- 
dure Q with the short-distance cutoff £, which translates 
into a derivative coupling in the limit of small £, 

4(0)V4(0<MO)V>£(£) eA^B^dAL- 

In the absence of derivatives, the four fermion term con- 
tributes T 2 to the inverse lifetime TjT 1 . The derivatives 
correspond to four powers of momenta close to the Fermi 
points in the scattering rate Tk.k'->p,p' oc (k—k') 2 (p—p') 2 , 
which translates into an additional factor of T 4 . Further- 
more, since at temperatures T <C 2j£ suppression of the 
conductance is entirely due to two-electron scattering, 
we expect the effective charge e* = S/(2\(Ib)\) obtained 
from a measurement of the shot noise S in the backscat- 
tering current (I B ) to be e* = 2e [H . 

For K < 1/4, the | A2 1 — > 00 fixed point is reached 
at T — and the system becomes insulating. The 
field 4>{x — 0, r) is pinned at the minima of the co- 
sine potential H 2 located at ±(2n + l)v / 7r/4 for A2 > 
and ±2ny / 7r/4 for A2 < 0, with n E Z. The conduc- 
tance is restored at finite temperatures by instanton pro- 
cesses corresponding to the tunneling between nearby 
minima separated by A0 = From the relation 

j e — edt4>/y/Tr between electric current j e and <j) field, the 
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charge pumped by a single instanton (Fig. \2b) is obtained 
as AQ inst = ;^=A</> = f • This fractionalized tunneling 
current can be understood as the Goldstone-Wilczek cur- 
rent 

[3 0, H3 for ID Dirac fermions with a mass term 
SC = g^>{mi + vy^m^)^ where g ~ A2, and the mass or- 
der parameters mi = cos 2y / 7r</>, m2 = sin 2y / 7r0 change 
sign during an instanton process with Acj> — i/n/2. The 
order is Ising-like because the 2P backscattering term 
explicitly breaks the spin U(l) symmetry of the helical 
liquid H = ^ J dx{Ka 2 z + j^p 2 ) down to Z 2 , where 
a z = p + — p- and p = p + + p_ are the spin and charge 
densities, and p± are the chiral densities for the two mem- 
bers of the Kramers pair. 

Fractionalization of the tunneling current is confirmed 
by a saddle-point evaluation of the path integral for 
large A2 in the dilute instanton gas approximation, which 
yields a Coulomb gas representation of the partition func- 
tion that can be mapped exactly to the boundary sine- 
Gordon theory 



S[6] = ^-y2\oj n \\0(^n)\ 2 +t f dr cosV^(t), 



where t is the instanton fugacity. The RG equation 
{}- — Tfc)t an d the conductance 



di 



2(l/4if-l) 



for 



for t follows as 

G{T) oc t{T) 2 is a power-law G oc (T/T|) 
T <C T£, K < 1/4. In contrast to the strong coupling 
regime in a usual Luttinger liquid where t corresponds 
to a single-particle hopping amplitude 22j, the unusual 
scaling dimension of the tunneling operator in the present 
case corresponds to half-charge tunneling. In particular, 
we calculate the shot noise in the strong coupling regime 
using the Keldysh approach and find S — 2e*|(J)| 
where (I) is the tunneling current and e* = e/2. 

Experimental realization. — We find that the experi- 
mental results of Ref. 7] are consistent with our the- 
oretical expressions for the weak coupling regime with a 
weak Luttinger parameter K ~ 1, but the small number 
of available data points does not allow for a reliable de- 
termination of the model parameters. The temperature 
dependence of the conductance being exponentially sen- 
sitive to K, our predictions can be best verified in QW 
with stronger interaction effects. Due to reduced screen- 
ing of the Coulomb interaction [24j , we expect to see a 
steeper decrease of conductance with decreasing temper- 
ature in HgTe samples with only a backgate. 

Because of lower Fermi velocities vf , we expect even 
stronger interaction effects to occur in InAs/GaSb/AlSb 
type-II QW which have been recently predicted to 
exhibit the QSH effect For QW widths Wi nAs = 

WGaSb = 10 nm in the inverted regime (26[, and consid- 
ering only screening from the front gate closest to the QW 
layer, from a k • p calculation of material parameters we 
obtain K ~ 0.2 < 1/4, making the insulating phase ob- 
servable at low temperatures. Although the backgate will 
cause additional screening, vp can be further decreased 



by adding a thin AlSb barrier layer between the InAs 
and GaSb QW layers. The Fermi velocity is controlled 
by the overlap between electron and hole subband wave- 
functions 3] which are localized in different QW layers in 
the type-II configuration (2r3 |. and an additional barrier 
layer will decrease this overlap. A lower vf also trans- 
lates into higher Kondo temperatures since vJ oc l/v F , 
where one power of vf comes from the matrix element of 
the localized impurity potential between edge states, and 
one power comes from the density of states v. Since Tj£ 
depends on vJ exponentially, we expect experimentally 
accessible Kondo temperatures in type-II QW. 
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